Exciton doublet in the Mott-Hubbard LiCuV04 insulator identified by spectral ellipsometry 
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Spectroscopic ellipsometry was used to study the dielectric function of LiCuV04, a compound comprised of 
chains of edge-sharing Cu04 plaquettes, in the spectral range 0.75 — 6.5 eV at temperatures 7 — 300 K. For 
photon polarization along the chains, the data reveal a weak but well-resolved two-peak structure centered at 
2.15 and 2.95 eV whose spectral weight is strongly enhanced upon cooling near the magnetic ordering temper- 
ature. We identify these features as an exciton doublet in the Mott-Hubbard gap that emerges as a consequence 
of the Coulomb interaction between electrons on nearest and next-nearest neighbor sites along the chains. Our 
results and methodology can be used to address the role of the long-range Coulomb repulsion for compounds 
with doped copper-oxide chains and planes. 



The elementary charge excitations observed by optical 
spectroscopy are a key source of information on electronic 
correlations in transition metal oxides. Because of their par- 
ticularly simple electronic structure (with a single hole in the 
d-electron shell) and their importance for high-temperature 
superconductivity and other electronic ordering phenomena, 
insulating copper oxides have long been recognized as model 
systems for charge excitations in the strong-correlation limit. 
In these systems, the electron-hole pairs created by optical ex- 
citation, — > form the Zhang-Rice-singlet 



(ZRS) state 111] with an oxygen-ligand hole in one CUO4 pla- 
quette and the full-shell configuration in a neighboring 
plaquette. In formal analogy to the single-band Hubbard 
model 10. S 01 J one usually refers to the continuum of ZRS 
excitations as the upper Hubbard band. Exciton formation of 
the {cPLh)i — df pairs has been considered as a possible 
low-energy elementary excitation in insulating copper-oxide 
planes or chains Qiilli. 

In analogy to semiconductors, excitons may form in 
narrow-band Mott-Hubbard (MH) insulators if the attraction 
between holes and electrons exceeds their kinetic energy. In 
copper oxides, the exciton size corresponds to the interatomic 
Cu-Cu distance, and hence a large binding energy of the or- 
der of an eV is expected, which is characteristic of a tightly 
bound, Frenkel-type exciton. In contrast to Frenkel excitons 
in semiconductors, however, electrons in the MH model can- 
not be localized on the same site because of the large on-site 
Coulomb interaction, U , and an intrinsic broadening of the 
exciton linewidth is expected as a result of exchange interac- 
tions 1 1 \\ . It has therefore been difficult to separate the in- 
trinsically broad and weak features due to MH excitons from 
the strong absorption continuum due to incoherent p-d tran- 
sitions il2lll3[l . Although peak structures arising from exci- 
tons have been identified in frequency-derivative reflectivity 
spectra 11311 . an accurate determination of the exciton spec- 
tral weight (SW) and a quantitative comparison with model 
calculations have thus far not been reported. 

In this Letter, we report a comprehensive ellipsometric 
study of charge excitations in LiCuV04, whose electronically 
active units are chains built up of edge-sharing CUO4 plaque- 
ttes. In this lattice geometry, the Cu-O-Cu bonds form ^ 90° 



angles, such that the nearest-neighbor inn) hopping amplitude 
is unusually weak and smaller than that between next-nearest 
neighbors (nnn). The resulting magnetic frustration and mul- 
tiferroic properties have recently drawn much attention to this 
family of compounds 1711 . In our optical exper- 

iment we observe a weak but well-resolved two-peak struc- 
ture for photon polarization along the chains. Using ellip- 
sometry, we accurately determine its temperature-dependent 
SW and show that it is controlled by the thermal variation of 
the nn and nnn spin correlation functions. The peaks can 
thus be identified as an exciton doublet arising from the long- 
range Coulomb interaction along the chains. By virtue of 
their exceptionally narrow electronic bandwidth, compounds 
with edge-sharing copper-oxide chains thus provide a highly 
favorable platform for the investigation of exciton formation 
and the interplay between spin and charge correlations in the 
cuprates. In the Cu02-planes of cuprate superconductors the 
band width is significantly larger, and hence excitons will not 
appear. However, a thorough understanding of the long-range 
Coulomb interaction is important for the quantitative descrip- 
tion of features such as charge density waves and stripes. 

LiCuV04 crystallizes in the orthorhombic space group 
Imma, with copper-oxide chains running along the 6-axis. Sin- 
gle crystals were grown from LiVOs flux, as described pre- 
viously iflS]. The magnetic susceptibility exhibits two char- 
acteristic maxima at T = 28 K, due to the establishment 
of magnetic coiTelations within the Cu-O chains, and at 2.4 
K, due to the formation of three-dimensional long-range or- 
der H [H I2III . For the optical measurements, the ab 
and ac surfaces were polished to optical grade, using a 0.25 
/im diamond suspension. The sample was mounted on the 
cold finger of a helium-flow cryostat with a base pressure of 
2 X 10~^ Torr at room temperature. Ellipsometric measure- 
ments were performed with a rotating-analyzer type Wool- 
lam VASE variable angle ellipsometer. Ellipsometry yields 
the anisotropic frequency-dependent complex dielectric ten- 
sor, = ei{bj) + ie2{'jj) = 1 + ! ATvaiuj) /uj, without the 
need for reference measurements or Kramers-Kronig trans- 
formations. A numerical regression procedure ll22ll was ap- 
plied to derive the principal components of the dielectric ten- 
sor from the ellipsometric data at different Euler angles. 
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Figure 1 shows the real and imaginary parts of the dielec- 
tric function of LiCuV04 at 300, 100, and 7 K. In order to 
separate contributions from the different bands to the optical 
response, and to explore the origin of the temperature depen- 
dence, we performed a classical dispersion analysis with a 
minimum set of Lorentzian oscillators by simultaneous fitting 
to ei{uj) and £2(1^^) 1I23I1 . One can distinguish the two most 
prominent optical bands located at 3.7 eV (3.55 eV) [3.6 eV] 
and 4.2 eV (4.4 eV) [4.4 eV] in the spectra along the a- (b-) 
[c~] axis. Clearly, the optical response is highly anisotropic. 
While the spectrum along the c-axis shows no significant tem- 
perature dependence, temperature dependent features are ap- 
parent along both a and b. The main temperature dependent 
features (marked by arrows in Fig. 1) are centered at 4.2 eV 
for polarization along a and at 2.95 eV along b. 

A detailed analysis of the spectra reveals that the origin of 
the temperature dependence for both polarizations is qualita- 
tively different. We first focus on polarization along a, per- 
pendicular to the chains. The temperature-difference spectra 
A(Ti(a>, T) and Aei{uj, T) (with consecutive intervals AT = 
50 K) displayed in Fig. 2(a) can be described mainly as a 
gradual narrowing and shift of the optical band at 4.2 eV with 
decreasing temperature, without any discernible change in its 
intensity. This behavior is typical for interband transitions and 
can be attributed to lattice anharmonicity. This conclusion is 
supported by the smooth temperature evolution of the af am- 
plitude at 4.2 eV and of e° at 3.8 and 4.45 eV displayed in 
Fig. 2(c). The only break in this gradual evolution is a subtle 
kink in the temperature dependence of cr° and £° at T ^ 30 
K [Fig. 2(c)], close to the maximum of the magnetic suscep- 
tibility. This may indicate some influence of spin correlations 
on the inter-chain charge transfer excitations. 

In contrast, the Aai{uj) difference spectra along the chains 
[Fig. 2(b)] show the emergence of a well-defined absorption 
peak at 2.95 eV at low temperatures, which is accompanied 
by an antiresonance feature with zero-crossing at the same 
energy in Ae5(t^). Figure 2(d) shows the temperature depen- 
dence of the crj amplitude near the center of this peak (upper 
trace) and the corresponding changes in measured at the 
off -resonant photon energies 2.65 and 3.45 eV (lower traces). 
These changes are Kramers-Kronig consistent and indicate a 
pronounced intensity enhancement of the band at 2.95 eV be- 
low 80 K, the temperature below which the magnetic suscep- 
tibility begins to deviate from the mean-field Curie-Weiss be- 
havior due to the appearance of short-range spin correlations 
along the copper-oxide chains | |21l| . 

The existence of positive and negative regions in the 
A(t"'''(w) spectra [Fig. 2(a) and 2(b)] indicate a redistribu- 
tion of optical SW. This can be quantified by integrating the 
optical conductivity in terms of the effective charge density 
AN,ff{u,T) = ^Jf^ Aai{to\T)dLo' , where m is the 
free electron mass and Ncu = 1.4 x lO^^cm"'^ is the density 
of Cu atoms. Figure 3(a) summarizes the low-temperature 
changes in SW. For polarization along the chains, the integral 
increases at low ui because of the band at 2.95 eV, but this 
gain is compensated by a SW loss within the spectral range 
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FIG. 1: Real S\(y') and imaginary £2(1^) parts of the dielectric func- 
tion of LiCuV04 measured at 7, 100, and 300 K for photon polariza- 
tions (a) E II a, (b) E |[ 6 and (c) E || c. Blue and red arrows mark the 
temperature dependent features discussed in the text. 



of the bands at 3.5 and 4.4 eV, such that the optical sum rule 
is satisfied. For polarization perpendicular to the chains, the 
integral of Act" iy) is consistent with a narrowing of the band 
at 4.2 eV upon cooling, and no SW redistribution between the 
3.7 and 4.2 eV bands is indicated. 

We now compare the salient observations to the results of 
model calculations, beginning with the overall anisotropy of 
the optical spectrum (Figs. 1 and 2). Since the main tran- 
sitions are expected to be due to Op^ Cud^ — > Op' Cud^" 
excitations, we have computed the optical spectrum of the 
previously proposed ll24ll pd-Hamiltonian, with on-site ener- 
gies ep — Ed ~ 3.5 eV and hopping parameters tpd = 1.4, 



pp 



0.55, tl 



0.5, t: 



2b 



3(b)]. While 



pp - ^-^.^pp - 0.4 eV [Fig. 
the overall agreement with the experimental data is reason- 
able, the mode at 2.95 eV is not described by the single- 
electron pd model, and the strong temperature dependence of 
its SW suggests that it is generated by many-body effects. Re- 
cent calculations of the optical conductivity of edge-sharing 
copper-oxide chains within a pd model including the on-site 
Coulomb repulsion 02511 have assigned the lowest-energy tran- 
sition along the chains to ZRS excitations. A strong tempera- 
ture dependence of the SW due to spin-correlation effects was 
also predicted, in qualitative agreement with our data. 

However, close inspection of the ellipsometric data indi- 
cates a more complex behavior that cannot be understood in 
terms of incoherent ZRS excitations. Fig. 4(b) shows that the 
2.95 eV band exhibits a double-peak structure with a satel- 
lite peak at 2.15 eV that is not visible in Figs.(l)-(3) because 
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FIG. 2; (a), (b): Temperature-difference spectra Aai{T,uj) = ai(Ti,Lu) - ai(T2,uj) (red) and Aei(r,tj) = £i{Ti,lu) - £i{T2,u!) (black) 
of LiCuV04 for polarizations (a) perpendicular to and (b) along the chains. Successive Aai (T, cj) and Aei (T, ui) spectra are shifted by 130 
D,~ cm" and 0.26 for clarity (2311 . The arrows mark the same energies as in Figs. 1(a) and 1(b). (c), (d): Temperature dependence of ai 
and si measured at 3.8, 4.21 and 4.45 eV for polarization perpendicular to the chains and at 2.65, 3.0 and 3.45 eV for polarization along the 
chains, as marked by vertical lines in (a) and (b), respectively. Cooling-down and warming-up curves are consistent and were averaged. 
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FIG. 3: (a) SW changes AN^ff{uj) ^ N^ff{Lu,50 K)-N^ff{uj,7 
K) in LiCuV04 for polarizations along (open squares) and perpen- 
dicular (open circles) to the chains. The arrows mark the same ener- 
gies as in Figs. 1(a) and 1(b). (b) £2 spectrum for polarizations along 
and perpendicular to the chains derived from experiment at 7 K and 
calculated from the pd-Hamiltonian as described in text. The optical 
bands above 5 eV listed in Table I were subtracted. 



of its low intensity II23I1 . The temperature dependence of 
the individual band SW was parameterized as A A^^^^ (T) = 
j^^^ASW'^^^T) in Fig.4(c), where ASW'^^'>{T) is deter- 
mined by the oscillator strength of the individual band. The 
SW of both bands increase upon cooling down to ^ 80 K in 
a parallel manner For lower temperatures, however, the in- 
tensity of the lower-energy band saturates, while that of the 
higher-energy band keeps increasing. 



We now demonstrate that the two bands can be quanti- 
tatively described as an exciton doublet in a model incor- 
porating long-range Coulomb interactions. We consider the 
one-band Hubbard model obtained by downfolding from the 
pJ-Hamiltonian \26] . In the case of edge-sharing chains, 
it is important to consider both nn (ti) and nnn {t2) hop- 
ping [Fig. 4(a)]. Besides the kinetic energy terms A',,; = 



the local Hubbard interac- 



tion U and the long-range Coulomb interaction Vi are in- 
cluded. The model thus reads 12711 : 



H 



(1) 



1=1 



Here Vi = V/l is, parameterized by the nn Coulomb interac- 
tion V . An optical excitation induced by the current operators 
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- h.c.) generates an empty site 
(ZRS) and a doubly occupied site (doublon in the upper Hub- 
bard band, UHB) one (di = h) or two (^2 = 26) Cu sites apart 
[Fig. 4(a)] iltl. As a result of the Coulomb attraction between 
the positive hole and the negative doublon, exciton states 
emerge below the UHB. For nn hopping only, the lowest- 
energ y ex citon dispersion is small, ^ — (2ti cos(/c/2))^/T^ 
our case, however, the dispersion is large and 
given hy Ei{k) = U — V ± 2t2 cos(fc) due to the t2 process, 
and the twofold degeneracy is lifted. The second exciton is 
centered at [/ — V/2. As seen in the inset of Fig. 4(c), the 
SW is concentrated in the exciton states, and only little op- 
tical weight is found in the UHB, despite its large density of 
states. While this representation gives a correct description 
of the energy scales, one would expect sharp excitons if one 
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FIG. 4: (a) Mott-Hubbard exciton formed by (fLh and 
states generated by nn (ti) or nnn{t2) hopping, (b) Temperature- 
difference spectra of £2(1^) near 2.9 eV for polarization along the 
chains, background corrected due to the temperature variation of the 
high-energy optical bands I23I1 . (c) Temperature dependence of SW 
of exciton bands at 2.15 and 2.95 eV. Inset: aiiuj) and density of 
states D{lo) of the UHB calculated with Eq. (1) at T = and pa- 
rameters U = 3.75, V = 1.6, ti = 0.08 and t2 = 0.1 eV. 



considered the charge sector only. We argue, however, that 
the coupling of the spin degrees of freedom (mainly due to 
i2-processes) implies a momentum average over the exciton 
band. This leads to a band with a width of ^ 4t2 ~ 0.4 eV 
for the lowest exciton, consistent with our experimental data. 

The SW for the individual exciton transitions can be ex- 
pressed in terms of the kinetic energy per bond and the hop- 
ping length di as AA^^^'^ = ~^df{Ki^i). The kinetic energy 
of the model given in Eq. (1) is equivalent to the superex- 
change energy of the corresponding Heisenberg model 112 811 
Hs = El Ji ■ S,+i - 1/4) - J[ ■ S,+i where 

J( ~ Atf /{U — Vi). Thus, the weights of the first and second 
excitons are directly related to the corresponding spin corre- 
lations: 



AN 



(0 



2m 
1^ 



■dfJi(S,-S,+i~l/A). 



(2) 



This relation explains why the second exciton is much 
stronger: (i) d2 = 2di and (ii) the nn spin-correlation func- 
tion {I = 1) is frustrated and small, because Ji — j[ is fer- 
romagnetic and small compared to the antiferromagnetic cou- 
pling constant J2. For / ^ 2, the correlation function is large 



(2) 

and negative, leading to large variation in AN^j:'^ as a function 
of temperature. These arguments also explain why the first ex- 
citon band experimentally observed at 2.15 eV is suppressed 
below 80 K, whereas the second band at 2.95 eV shows a steep 
increase [Fig. 4(c)] following the spin-correlation effects. A 
numerical computation of Eq. (2) using exchange parame- 
ters determined by neutron scattering 11911 estimates within a 

(2) 

factor of two the experimentally determined AN^j:'j:{75 K— 7 
K) ~ 0.013. Within the model the exciton mode draws its SW 
from the UHB, which may be superimposed on the charge- 
transfer excitations at 3.5 and 4.4 eV [Fig. 3(a)]. 

In summary, our optical measurements have yielded clear 
evidence for exciton formation and interesting insights into 
the relationship between the charge dynamics and frustrated 
magnetism in a Mott-Hubbard insulator with edge-sharing 
copper-oxide chains. The results and methodology estab- 
lished here are a good starting point for further investigations 
of compounds with doped copper-oxide chains and planes. 
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Anisotropic dielectric response of LiCuV04: overall features and dispersion analysis 
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FIG. SI: Real ei{u>) and imaginary e2(<^) parts of tlie dielectric function of LiCuV04 measured at 7K for photon polarizations (a) E || a, (b) 
E II & and (c) E || c and represented by the individual contributions of separate Lorentzian bands. The parameters are listed in Table SI. 



In order to separate contributions from the different bands to the optical response of LiCuV04, we per- 
formed a classical dispersion analysis by simultaneous fitting to £i{uj) and e2{uj). The decomposition of the 
optical response in terms of a sum of Lorentzians is a phenomenological fitting procedure to quantitatively 
parameterize the observed temperature-driven optical anomalies and spectral weight transfer. 

Moreover, the agreement between the simulated and measured spectra confirms that our data for ei (lu) 
and €2{uj) (or ai{uj)) determined independently by spectroscopic ellipsometry are Kramers-Kronig consis- 
tent. 



TABLE S 1 : Parameters of Lorentz oscillators resulting from a dispersion analysis of complex dielectric response in a-axis (&-axis) [c-axis] 
polarization in LiCuV04 measured at T = 7 K. eoo=l-75 (eoo=l-69) [eoo = l-55]. 

LOjieV) Sj (eV^) FjCeV) 

3.31 (2.93) 1.86 (0.52) 0.62 (0.47) 

3.70 (3.55) [3.59] 9.30 (5.92) [6.83] 0.67 (0.93) [1.11] 

4.20(4.41) [4.38] 20.6 (7.19) [15.7] 1.07 (0.81) [2.03] 

5.94 (5.47) [5.58] 21.5 (16.8) [2.80] 1.62 (1.76) [1.11] 

6.86 (7.00) [7.00] 24.0 (66.2) [59.3] 1.00 (1.90) [2.79] 
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A minimum set of Lorentzian oscillators, with one high-energy oscillator beyond the investigated spectral 
range, was introduced to represent a dielectric function in the form e{uj) = £00 + J2j uj'^~oj^^~iujr ' where 
ujj, Tj, and Sj are the peak energy, width, and oscillator strength of the jth oscillator, and £00 is the core 
contribution from the dielectric function. The parameters determined by simultaneous fitting to £i{uj) and 
£2(1^) measured at 7 K are listed in Table SI. 

• Temperature-difference spectra 

Figure S2 shows the temperature-difference spectra Acri(T, uj) = ai{Ti,Lu)—(7i{T2, to) and Aei(T, u) = 
ei{Ti,u) — Ei(T2, uj) with consecutive intervals AT = 50 K and represents an expanded version of Figs. 
2(a) and 2(b) of the manuscript. 




FIG. S2: (a), (b): Temperature-difference spectra Acri (T, cu) = ai {Ti,lu) — ai(T2,Lj) (red) and Aei{T, u) = ei(Ti, a;) — ei{T2,u)) (black) 
of LiCuV04 for polarizations (a) perpendicular to and (b) along the chains. Successive Acti(T, oj) and Aei{T,uj) spectra are shifted by 
130 Sl^^cm"^ and 0.26 for clarity. The arrows mark the same energies (4.2 eV and 2.95 eV, respectively) as in Figs. 1(a) and 1(b) of the 
manuscript. 



• Optical response along the chains: Double-peak structure near 2.9 eV 

Figure S3 shows that the emergence of the satellite absorption peak at 2.15 eV at low T is evident from 
the obvious changes in both the imaginary, e2{to), and real, £1(0;), parts of the complex dielectric function 
presented in Fig. 1 (b) of the manuscript. The arrow in the lower panel marks, in particular, the photon energy 
at which the positive changes in 62{uj) coincides with crossing in ei{u), as temperature decreases from 300 
K to 100 K. These features are clearly resolved in the temperature-difference spectra shown in Fig.S4, as 
a "bump" in Ae2{uj) and a "wiggle" in Aei{uj), and can be quantitatively parameterized by the Lorentzian 
fit. From this fit, the temperature dependence of the individual band spectral weight is determined by the 
oscillator strength of this band, ASW'^'-^ (T) = Si{T) [eV'^]. Alternatively, one can integrate the real part of 
the optical conductivity within the individual band, A S'IV('^(T) = l/An J*^'^ a;Ae2(i^)rfi^, after subtraction 
of the electronic background due to the temperature variation of the high-energy optical bands. The gradual 
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FIG. S3: Real ei(cj) and imaginary £2(1^-') parts of tlie dielectric func- 
tion of LiCuV04 measured at 7, 100, and 300 K in 6— axis polariza- 
tion. Upper panel is the same as in Fig. 1(b) of the manuscript. Lower 
panels show an enlarged view of the upper panel. 



FIG. S4: Temperature-difference spectra of imaginary e2(<^) (upper 
panel) and real ei{uj) (lower panel) parts of the dielectric function 
of LiCuV04 for polarization along the chains. Vertical dashed lines 
mark photon energies 2. 15 eV and 2.95 eV at which resonant behavior 
in Ae2{u!) coincides with zero crossing in A£i(a;). 



changes of the electronic background with temperature are well controlled by the dispersion analysis of the 
whole £1(1^) and e2{uj) spectra performed at every temperature. Figure 4(b) of the manuscript shows the 
background corrected temperature differences spectra of e2{uj) near 2.9 eV. Both the above procedures 
give the same temperature dependence of the individual band spectral weight, the uncertainty is within the 
symbol size in Fig.4(c) of the manuscript. 



